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The behaviour of a relativistic electron-positron plasma in the pulsar magne- 
tosphere is investigated. The equation of the motion of the magnetospheric 
plasma is discussed, from which it follows that, if the plasma particle radial 
velocity > c/\/2 (where c is the speed of light), the centrifugal acceleration 
^ ■ changes its sign and the particle braking begins. The stability of the magneto- 

spheric plasma with respect to the radially oriented potential perturbations is 
discussed and the possibility of the electrostatic field generation in this plasma 
^ ■ along the pulsar magnetic field lines is shown. 

^ ■ PACS numbers: 97.60. Gb, 94.30.-d, 94.30. Cm, 94.30.Kq 
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I. INTRODUCTION 



The investigation of the pulsar magnetosphere is of great interest (see, 

for 

example [1-9]). For the determination of the magnetosphere structure it is 
very important to study the physical processes in plasma of magnetosphere. 
After the pioneer papers [1,2] it is assumed that, because of the pulsar coro- 
tation with its magnetic field, the electric field is generated, which has the 
nonzero component along the magnetic field. The electric field ejects parti- 
cles from the pulsar surface and accelerates them up to relativistic velocities. 
The particles moving along curved magnetic field lines radiate 7-quanta and, 
when their energy exceeds electron's doubled rest energy 2mc^ {e^>2mc^), 
7-quantum decays into an electron-positron pair. This pair is also acceler- 
ated in the electric field and 7-quanta appear again, which again decay into 
an electron-positron pairs, etc. This cascade process causes a fast filling of 
the magnetosphere with the relativistic electron-positron plasma, which, in its 
turn, screens the electric field generated by the pulsar rotation. In the previous 
papers (see, for example [4,5]) the strongly turbulized electron-positron plas- 
mas were investigated, but the process of the turbulization was not discussed. 
We investigate the possible ways of the plasma turbulization, so we discuss the 
processes in the approximation of the weak turbulence. 

In the first section of the present paper the equation of the motion for the 
magnetospheric plasma in the zeroth approximaton of the weak turbulence 
is discussed, from which it follows that, if the plasma particle radial velocity 
Vr > (where c is the speed of light), the centrifugal acceleration changes 

its sign and the particle braking begins. 

In the second section the equations for the perturbed quantities (in the first 
approximation of the weak turbulence) are presented. Also the stability of the 
magnetospheric plasma with respect to the radially oriented potential pertur- 
bations is discussed and the possibility of the electrostatic field generation in 
this plasma along the pulsar magnetic field lines is shown. 
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II.THE EQUATION OF THE MOTION OF THE 
MAGNETOSPHERIC PLASMA 



The equation of the motion of the magnetospheric plasma particles was 
discussed in paper [10]. We use the perpendicular rotator model of the pulsar 
magnetosphere and treat only the polar cap. The magnetic field lines are con- 
sidered as the radial straight lines located in the plane, which is perpendicular 
to the pulsar rotation axis. This assumption is justified, because we discuss 
the processes in the magnetospheric layer, the thickness of which is much less 
than the curvature radius of the magnetic field lines. The magnetospheric 
plasma particles move along the pulsar magnetic field lines and also corotate 
with them, because the field lines are frozen in the plasma. The electric field, 
generated by the pulsar rotation together with its magnetic field, is screened 
by the magnetospheric plasma. 

It is convenient to begin the discussion of the plasma particle motion in the 
noninertial frame of a rotating magnetic field line, which is described by the 
metric: 



dS^^-{l-Q^r^)dt^+dr^, (1) 

where Q is the pulsar rotation frequency. Here and below we use so called 
"geometric units": c — G—1. 

According to the Einstein principle of equivalence, we can not tell gravi- 
tation from noninertiality. Thus, for the description of the particle motion in 
the pulsar magnetosphere the "3-1-1" formalism can be used. This formalism 
is described in [11]. According to the "3-1-1" formalism, the equation of the 
motion for the particle with the mass m and charge e has the following form 
[11]: 

-%+{VV)p^-^—+-{E+[VB]), (2) 

a ot am 



where a is the so called " lapse function" and in our case a = — fi^r^. Here 
and below we use the dimensionless momentum p {pis changed by p/m). We 
can rewrite the equation (2) for the quantities defined in the rest inertial frame: 

^+{VV)p=-^aVa+-{E + [VB]). (3) 
ot m 



3 



Now let us discuss the motion of the plasma particles in the zeroth appro- 
ximation of the weak turbulence. In the limits of this approximation the 
quantities, which are located in the equation of motion, can be presented as: 



E^Eo+Ei, B^Bo+Bi, p^p'o+p'i, (4) 

— # — # — * — * 

where £"0; Bq and po are the basic terms and Ei, Bi and pi are the perturba- 
tions in the first approximation of the expansion over the parameter of the weak 
turbulence. The small parameter in the approximation of weak turbulence for 
the electron-positron plasma is: 

(5) 

mn'f 

Prom the equation of the motion in the zeroth approximation for the radial 
acceleration one can obtain: 

^+(V^VW=-7o^^V. (6) 

Let us introduce new variable I, which is connected with coordinate r by 
the following relation: 

r^l+ f v{l,T)dT. (7) 

J to 

It is clear, that / is a coordinate derived at the moment to- Thus, doing the 
transformation to the Lagrangian variables and using equation (6), we receive 
(see also [12]): 



(8) 



Equation (8) can be solved exactly. Using Jacobian functions, the solution 
can be presented in the form [12]: 

where Sn and dn are the Jacobian elliptical sine and modulus respectively 
[13], Voi is the particle initial velocity. From equation (8) it follows that, if the 
radial velocity Vr > 1/a/2, the acceleration changes its sign and the particle is 
not accelerated, but braked (see also [12]). 
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In the case Voj— >1, using the asymptotic expression for the Jacobian func- 
tions one can find [12]: 



l{t)^-^sinnt. 



(10) 



For the radial velocity we will obtain: 



Vor=Voj^cosilt, 



(11) 



from which it follows that 



Vo''^{Vorf+{Vo^f^const (12) 
(because of the corotation, Vq^ = ilr), i.e. no energy is expending on the 



particle braking along the field line, the energy transforms from radial to the 
transversal one. 

III. THE EQUATIONS FOR THE PERTURBED QUANTITIES 

As it was shown above, the relativistic plasma particles are braked in the 
pulsar magnetosphere, if their radial velocity > 1/ y/2. It is very interesting 
to discuss the stability of such a plasma with respect to the radial perturba- 
tions. In particular, we discuss the potential perturbations oriented along the 
magnetic field lines. The initial stage of the perturbation development can be 
described by the equation, which is easy to obtain from (3) by substituting in 
it the expansion (4). For the first order terms one can obtain: 



where ni—nip—riie is the perturbed concentration — the difference of positron 
and electron perturbed concentrations. 

Let us complete the equations (13) and (14) with the continuity equation 
in the first approximation of the weak turbulence and thus, lock the set of the 
equations: 




divEi=4:7Ten 



(14) 



5 



dm -> -> 

— — \-divnQVi-\-divniVQ—Q. (13) 
ot 

Before solving the set of equations (13)-(15), let us introduce the charac- 
teristic spatial parameter of the system and consider the case, when 

[^M.]'' » A, (16) 

where A is the length of the perturbations. Using the value of Vqj. obtained 
from equation (12), equation (16) could be presented in the form: 

» A, (17) 

where Rlc — C/Q,\s the radius of the Light Cyhnder. It follows from equation 
(17) that, if we investigate the processes in the region with the scale smaller 
than Rlc one, we can consider the medium as homogeneous one and make the 
Fourier transformation. Vq - is solved in our description in Lagrangian vari- 
ables, while the set of equations (13)-(15) is written in the Eulerian variables. 
So, let us express Vqj. by the variables r and t (see (7)). Doing the spatial 
Fourier transformation using approximation (16) for the perturbed quantities 
Pi, rii and Ei, we can obtain equation, which corresponds to the set of the 
eqs. (13)-(15): 



— +lkrVor 

at 



PirSin2D,t —pirSin'^2D,t, (18) 

7o 



where Up is the plasma frequency. 

In order to solve the equation (18), it is convenient to introduce an expres- 
sion: 

krR{t)^ J krVor{t)dt. (19) 
Then the equation (18) can be rewritten in the following form: 

-^{exp{ikrR)pir {r,t))^ — {exp{ikrR)pir{'r,t)sin2nt) - 
(J ^ 

—exp{ikrR)pir{r,t)sin'^nt. (20) 
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Let us expand the function 
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where 



exp{±ikrR)—exp{±iasinD,t), 



(21) 



a~, (22) 



in the row over the Bessel functions: 



+00 



exp{±iasinflt)= ^ In{a)exp{±inflt). (23) 

n=— 00 

Here /„(a) is a Bessel function. Doing the expansion (23) in equation (20) and 
using the Fourier transformation in time, we obtain: 

^2 +00 

-^Pi^(f,a;)= Is{a)In{a){'^-n^fpirir,uJ+{s-n)n)- 
^7o 



s,n=— 00 



^ Is{a)In{a){u;—nO,)pir{f,u;+{s—n+2)Q)+ 
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s,n=— 00 



Is{a)Inia){uj-nfl)pir{f,uj+{s-n-2)Q)- 
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s,n=— 00 

+^Pi,(f,(a;-21]))+^Pi,(f,(a;+20)), (24) 
470 470 

where cu is the frequency of the generated perturbations. We will consider the 
case of very low frequencies ux^ft, which seems to us to be the most interes- 
ting. We will assume that the particle has no time to change its momentum 
significantly, which fully corresponds to the linear approximation of the weak 
turbulence and describes the initial stage of the process. In this case all har- 
monics except the zeroth pi^(r,a;) are small. So keeping only zeroth harmonics 
in eq.(18), which corresponds to the averaging by the fast oscillations and 
changing n— >s, we obtain the dispersion relation: 

2 +00 QTr 2 p +00 



s=— 00 



270 s^oo 4 
+00 +00 

-D, Y sIs{a)Is-2{a)-{ui-2Q) Y Is{a)Is+2{a)+^l Y sls{a)ls+2{a) 

s=—oo s=— 00 s=—oo 

It is easy to find that: 
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+ 00 +00 



s=— oo s=— oo s=— oo 

+00 +00 2 

5: /,2(a)=l and ^ «'/.'(a) = y. 

s=— 00 s=— 00 

So from the (25) we can obtain: 

" -2^ 

We know that Eir^exp{~iujt), so one can conclude that, when the second 
term in the right hand side of (26) is larger than the first term the aperiodic 
instability is being developed in the pulsar magnetosphere, i.e. the field Ei is 
increasing exponentially along the magnetic field lines. 

The condition of the aperiodic instability development can be written in 
the following form: 



< (27) 



Up 



where Ir is the charge separation scale in the magnetospheric plasma. For the 
typical parameters of the pulsar magnetosphere the charge separation scale at 
the light cylinder (the light cylinder is the surface, on which the azimuthal 
velocity equals to the speed of light — Qr — c) is of the order lO^sm. As 
for the parameter (Voj^/7o)/a;p, it differs from the relativistic generalization of 
the Debye radius Id [14] {Id — (yTy/lr) /'^p, where Vt and 7^ are the particle 
thermal velocity and Lorentz- factor respectively). On the other hand, it is 
selfevident that the following condition must be fulfilled: 

lr'>lD, (28) 

i.e. the charge separation scale Ir must be much larger than the Debye radius 
Id- Prom the conditions (27) and (28) one can conclude that the aperiodic 
instability will take place if: 

7o>7oT- (29) 
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IV. CONCLUSION 



At the end let us discuss the possible results of the instability. We can 
see that the plasma motion along the magnetic field lines and at the same 
time rotation together with them (i.e. corotation) cause the generation of 
the aperiodically increasing electrostatic field under the condition (27). On 
the other hand, it is selfevident that the corotation can not take place on 
the arbitrary distances from the pulsar surface, because on some distance the 
azimuthal velocity will reach the speed of light = Qr = c. So, the corotation 
must be removed. The instability, which was discussed above, can contribute 
to the process of the corotation removing, in particular, the increasing electric 
field will cause the additional braking of the particles of one sort and the decrea- 
sing of braking of the other sort. This fact will evidently cause the motion of 
the electrons and the positrons with respect to each other, i.e. the increasing 
current j will appeare. So, according to the Maxwell equation 47rj=rotB, 
the magnetic field will be generated. The current will be directed along the 
pulsar magnetic field lines, therefore, the generated magnetic field will have 
an azimuthal component B^p. The particles move along the field lines, so, the 
corotation will be removed. The electric field, i.e. the current j, will increase 
untill the corotation law removal. 
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